Exponential localization of Wannier functions in insulators 
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The exponential locahzation of Wannier functions in two or three dimensions is proven for all 
insulators that display time-reversal symmetry, settling a long-standing conjecture. Our proof relies 
on the equivalence between the existence of analytic quasi-Bloch functions and the nullity of the 
Chern numbers (or of the Hall current) for the system under consideration. The same equivalence 
implies that Chern insulators cannot display exponentially localized Wannier functions. An explicit 
condition for the reality of the Wannier functions is identified. 
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Wannier functions play a fundamental role in the de- 
scription of the electronic properties of solids [l| . They 
allow for an intuitive interpretation of the bonding prop- 
erties of solids [2], they are at the center of the mod- 
ern theory of polarization and they form a very ef- 
ficient basis for order-TV calculations or the construction 
of model Hamiltonians 0]. 

For onc-dimensional systems, Kohn |5:] proved that 
Wannier functions are exponentially localized. This 
property has many desirable consequences, such as the 
existence of moments (r") for all n, the exponential con- 
vergence of numerical calculations [6y , and the possibility 
to use Wannier functions for the description of surfaces 
0. In two and three dimensions, the existence of expo- 
nentially localized Wannier functions is one of the few un- 
solved problems of one-particle condensed-matter physics 
0. In the absence of a proof for exponential decay, this 
property is commonly used as working hypothesis [9^, 
checked with numerical calculations ^lOj , or simply taken 
for granted 11, 1^ . 

In this paper, we demonstrate that exponentially lo- 
calized Wannier functions exist for insulators in two and 
three dimensions, and show how localization is related to 
the Berry connection for the set of bands under consid- 
eration ^2] and to the corresponding Chern number(s). 
Chern numbers have been playing a rapidly increasing 
role in solid-state physics, from metal- insulator transi- 
tions [isj . to quantized transport jl4i] . transition metal 
nanomagnets [3], quantum Hall effect and its spin 
analogue fl7'|. Our central result is that Wannier func- 
tions with exponential decay can be constructed if and 
only if all the Chern numbers are zero. This implies that 
if the system is symmetric for time-reversal, then the 
Wannier functions are exponentially localized. 

We consider a crystal with a Bravais lattice F and a 
unit cell C. The reciprocal lattice is denoted by F* and 
the reciprocal unit cell by SI. A function /(r) is called 



periodic if /(r + R) = /(r) for any vector R of F, while a 
function /(k) is called periodic in the reciprocal space if 
f(k+G) = /(k) for any vector G of F*. The dynamics of 
the electrons in the crystal is described by a Hamiltonian 
H — —A + V{r) (in Ry), where the potential V is real 
and periodic. According to Bloch's theorem, the eigen- 
functions of H can be written as ■!/'nk(r) = e*'' '"u„k(r), 
where the Bloch states Urak(r) are periodic eigenstates of 
the Hamiltonian H{'k) = (— iV -I- k)^ -|- ^(r) with energy 
e„(k) and satisfy the boundary conditions 



Unk+G 



(r) 



-iG r 



(r). 



(1) 



Wannier functions are defined as 



z«„(r-R) = ±^dke'''-(--^)u„k(r), (2) 

where \ is the volume of Q. The localization properties 
of the Wannier functions are related to the regularity of 
as a function of k. In a nutshell, the more regular 
the states, the more localized the Wannier functions [^, 
18 . T^ . Exponential decay is obtained if and only if the 
functions are analytic [3, 0] ■ 

In the simplest procedure, the energies e„(k) and 
the Bloch functions u„k are determined as a function 
of k by ordering the eigenvalues by increasing energies 
< e2(k) < .... Although this procedure is stan- 
dard in band-structure calculations, it gives Bloch states 
Unk that are not regular in k since the phase of u„k is 
random. Moreover, at a crossing point the energy e„(k) 
can have a kink. In one dimension, Kohn Q showed that 
it is possible to define Bloch states that are analytic func- 
tions of k. In two and three dimensions this is generally 
not possible 18| because band crossings lead to Wannier 
functions that can decay as [2l|. 

Blount [24I noticed that the decay properties of Wan- 
nier functions can be improved by considering a set 
of eigenstates, called a composite band [23| . which are 




FIG. 1: Band structure of silicon. The contour 7 encloses the 
energies of the composite band. 



seaparated by a gap from all others. More precisely, 
if ei(k), . . . , eA/(k) are the eigenstates of the compos- 
ite band, we assume that there is an a > such that 
|e„(k) — em(k)| > a for any k in 17, 1 < n < A/ and 
m > M. In insulators it is always possible to define 
a composite band since gaps exist between valence and 
conduction bands and between valence bands and core 
states. This is not necessarily the case in metals. An 
example of a composite band is given in Fig[T]for silicon. 
Exponentially localized Wannier functions are obtained 
if we can find M analytic functions D„k (called quasi- 
Bloch states) that span the same vector space as the M 
Bloch states of the composite band. The problem of the 
existence of analytic quasi-Bloch states was solved for a 
sing le isolated band (i.e. M = 1) in 3D by des Cloizeaux 
[iSj and Nenciu In this work we determine when ana- 
lytic quasi-Bloch functions exist for any composite band 
in two and three dimensions. 

The quasi-Bloch functions w„k can be expressed as 
I'nk = Wmk(r)?7mn(k), whcrc C/(k) is a unitary ma- 
trix, defining the Wannier functions 



(r)[/„„(k). (3) 



The improved localization of these Wannier functions was 
observed on many systems 0, [IJ] . 

To study the analytic properties of the Bloch functions, 
we consider a complex vector k = k' + ik". Then, -ff (k) 
is not Hermitian but there still are eigenvalues en(k) 
and eigenstates w„k such that i?(k)it„k = e„(k)u„k. 
The branch points of en(k) determine the points of non- 
analyticity of the function u„k @, . Even very simple 
crystals, such as silicon, exhibit such branch points [25 1. 

If the potential ^(r) is square integrable, it can be 
shown that i?(k) is analytic [26]. The reality of V{r) im- 
plies that iJt(k) = H(k*) and H*(k) iJ(-k*). From 
a practical point of view, the condition of square integra- 
bility of V{r) encompasses potentials with Coulomb sin- 
gularities and the potentials used in LDA and GGA cal- 




FIG. 2: Mobius strip. The solid circle is the base B and 
the thin straight lines are the fibres F{x). The dashed line 
represents the constant but not antiperiodic basis e. The 
arrows represent the continuous and antiperiodic basis vectors 
v(a;). It is clear that v(a::) has to be zero for some x. 



culations. The calculation of band structures in the com- 
plex plane is available in some standard band-structure 
packages [27 1. 

Several authors [1, [3, [111 noted that the obstacles to 
the existence of quasi-Bloch states are topological. To 
illustrate this, we introduce a simple mathematical con- 
cept, that of a fibre. It is quite common in physics to 
consider a space (let us call it the base) and to asso- 
ciate a vector space (called the fibre) to each point of the 
base. For example, in magnetostatics, the base B is the 
space M'^ and to each point r of the base we associate 
a three-dimensional vector space (the fibre F{r) = M'^). 
The vector potential A(r) is a vector of the space F{r). 
In the Born-Oppenheimer approximation, the base B is 
the set of possible positions of the nuclei Ri , . . . , Rat and 
the fiber _F(Ri, . . . , R^v) corresponding to a nuclear con- 
figuration is the vector space generated by the solutions 
of the Schrodinger equation with clamped nuclei. For a 
composite band, the base is the Brillouin zone f2 and the 
fibre F(k) is the Af -dimensional vector space generated 
by the Bloch states u„k- The quasi-Bloch states Vnk can 
now be defined precisely as a basis of i^(k) such that each 
Vnk is analytic in k, periodic in the reciprocal space and 
satisfies the boundary conditions (jlj. 

To show that topology might forbid the existence of 
quasi-Bloch states, we consider the simple example where 
the base B is a circle represented by the angular variable 
X € [0, 2tt] and the fibre is F{x) = R. We denote by \r{x) 
a basis of F{x). Topology intervenes when we determine 
how the fiber at a; = is related to the fiber at a: = 27r: 
we can use periodic boundary conditions v(27r) = v(0) 
(cylinder) or antiperiodic boundary conditions v(27r) — 
— v(0) (Mobius strip, illustrated in Fig. 2). In the case 
of a cylinder, we can choose a basis v(a;) = e (a non-zero 
vector independent of x). This basis is obviously regular 
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and periodic in x. In the case of a Mobius strip, let us 
assume that a regular and antiperiodic basis v(a;) of 
exists. Let y{x) be the component of v(a;) with respect 
to the periodic basis e, namely v(x) — y{x)e. Clearly 
y{x) has to be regular and antiperiodic (2/(0) = — 2/(27r)). 
By the intermediate value theorem there needs to be a 
point xq between and 27r such that y(xo) = 0, as shown 
in Fig. [2] Thus at xq, v(xo) — 0. But the null vector 
cannot be a basis for F{xo). Therefore, the topology of 
the Mobius strip implies that no regular basis of F{x) 
can exist. 

The topological obstruction to the existence of quasi- 
Bloch states in dimensions one, two and three has been 
recently studied j29i] . It was discovered that quasi-Bloch 
states exist if and only if all the Chern numbers of the 
system are zero. Such result is of central importance to 
our paper, since it shifts the focus of our analysis on the 
determination of the nullity of Chern numbers in dimen- 
sions two and three. 

For a given k, P(k) = Y^tLi |unk)('Unk*| defines a pro- 
jector 221 . The Riesz formula for the projector is 



P(k) 



1 

27ri 



dz 



z - H{k) 



where 7 is a contour enclosing all the eigenvalues e„k of 
the composite band (see figll]). P(k) is analytic in k on 
a strip f)c = {k = k' ik" ,k' e \k'J\ < AJ = 1, 2, 3} 
with A > even if the states u„k a-re not analytic in k 
[ist . The value of A is related to the band gap An 
example of contour 7 for Si is given in Fig. [T] 

To calculate the Chern numbers, we introduce the 
Berry connection corresponding to the basis u„k of F(k), 
defined in dimension d as the d-dimensional vector of ma- 
trix functions A(k) = (A^ (k), . . . , A'^(k)) with matrix 
elements 

A„„(k) = / <„k(r)VkM„k(r)dr. 
Jc 

The trace of the curvature of this connection is 0, [Isj 

[A\A^]) (4) 



/ dA^ dA' 



= Tr P(k) 



dP{k) dP{k) 



dki ' dkj 



where tr is the matrix trace and Tr the operator trace. 

In two dimensions, this curvature leads to the unique 
Chern number 

Ci = ^ I B^^{k)dk^dk2. (5) 

In three dimensions, we have three Chern- like numbers 

Ce ^ f B'^k)dhAdk,, (6) 



where £ = 1,2, 3. The domain of integration Qi is a torus 
defined as the set of points (fci, fc2, fcs) of such that ki = 
0. The reality of V{r) implies the time-reversal symmetry 
P(k) = (P(-k))*. Therefore, B'^ {-k) = -B'Hk), the 
Chern numbers are zero and quasi-Bloch functions exist. 

With this result in hand, we can repeat the reasoning 
of des Cloizeaux (section III.B of Ref. Ts^, to show that 
lim|i^|^oo e^l^l'u;„(r - R) = for any b < A. 

Our proof was given for the case of "spinless" electrons. 
To take spin, spin-orbit and all relativistic corrections 
into account, we consider a crystal described by the Dirac 
Hamiltonian ^ H{k) = ~ica- {V + ik) + f3c^ /2 + V{r). 
The Bloch functions are Dirac spinors u"y^{r) and the 
corresponding Berry connection is 

A-Jjk) = |dr«,(r))*Vku:^k«- 

The curvature is defined by equation Q , where the trace 
is over the n and a indices, and the Chern numbers by 
eqs.© and Again, quasi-Bloch functions exist for 
the Dirac Hamiltonian if and only if all Chern numbers 
are zero (in particular, if the potential V is square in- 
tegrable and time-reversal symmetric). In that case, the 
relativistic Wannier functions are exponentially localized. 

Given the existence of exponentially decaying Wannier 
functions u„k, we need an algorithm to determine the uni- 
tary matrix U{k) such that u„k(i") = J2m t^mn(k)Mmk(r)- 
We sketch now a possible approach. Being analytic, the 
quasi-Bloch functions satisfy the Cauchy-Riemann equa- 
tion Bvnk — 0, where B = {Bx, By, Bz) and, if k = k'-l-ik", 
Bj — [\/2){dk'. + idk")- If we denote by u(k) and v(k) 
the vectors with components Unk and w„k, the Cauchy- 
Riemann equation gives us Bv = {Bu)U + u{BU) = 0. 
Multiplying by the conjugate of u, we find BU = XU, 
where Xmn{k) = — (M^_kl5w„k)- Note that, on regions 
where w„k is analytic, Xmn{k) = 0. The equation BU = 
XU should be solved in the space of matrices U (k) peri- 
odic in the reciprocal space such that U{k)'' = U{k*)~^. 
This condition ensures the orthogonality of the quasi- 
Bloch functions on the real axis. Explicit integral expres- 
sions are available to solve the Cauchy-Riemann equation 
Bf = g [32]. They turn BU = XU into an integral equa- 
tion that can be solved numerically. 

Note that, if t/ is a particular solution of BU — XU, 
the general solution can be written UA, where A is any 
analytic matrix periodic in the reciprocal space and sat- 
isfying A(k)'l' — yl(k*)~^. For instance, A(k) can be any 
matrix of the form A{k) ~ exp ( an.e*'''^) , where 
the sum is over a finite number of sites of F and the 
matrices or satisfy = — a_R. This very large set of 
solutions corresponds to the fact that exponential decay 
is a long-range property: Any finite linear combination 
of exponentially decaying Wannier functions centered on 
various sites is still exponentially decaying. Therefore, 
it is still necessary to optimize localization around the 
centers of the Wannier functions 2] by properly choos- 
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ing ^(k). Such an approach provides Wannier functions 
that are locahzed on the short and long range. 

It is convenient to determine under which condition 
the Wannier functions are real. In eq. we take the 
complex conjugate and change the variables k — > — k 



w, 



;(r) 



-y 



dke^'^-<„_k(r)(C/„„(-k))' 



The reality of V{v) implies that we can choose u^^kC'') — 
Wmk(r). Thus, the Wannier functions are real if 
(^^mn(— k))* — L/m„(k) for real k. This property is sat- 
isfied if U satisfies L/(-k*) = [/(k)* for complex k. 

Wannier functions are often obtained by minimizing a 
functional Q.{JJ) 0. The time-reversal symmetry / acts 
as (/C/)m„(k) — ([/„!„(— k))* for real k. If the func- 
tional r2 satisfies the symmetry il{IU) = n{U) and if 
has a unique minimum (up to a possible overall phase), 
then U — lU and the Wannier functions are real. The 
spread functional f2 defined by Marzari and Vanderbilt 
[2| satisfies the symmetry n{IU) = Q{U). Thus, when 
the spread functional has a unique minimum f/, the cor- 
responding Wannier functions are real, proving for this 
case the conjecture of Ref. S 

In this paper, we assumed that the system is time- 
reversal symmetric. If this is not the case (e.g. for the 
Haldane Hamiltonian [111), the present approach implies 
that exponentially localized Wannier functions exist in 
regions of the parameter space where the Chern num- 
bers are zero or, equivalently, when the Hall current is 
zero. Thus, the vanishing of the Hall current is a measure 
of the exponential localization of the Wannier functions. 
This confirms rigorously Thouless' observation 3J|. As 
a corollary, we deduce that no exponentially localized 
Wannier functions exist for Chern insulators (i.e. insu- 
lators with non-zero Chern numbers (s^l) and that time- 
reversal symmetric systems cannot be Chern insulators 
(as noticed by Haldane [s^l for the case M = 1). 

In conclusion, in this work we demonstrated that Wan- 
nier functions are exponentially localized for insulators 
that satisfy time-reversal symmetry, and we showed that 
the vanishing of the Chern numbers is equivalent to the 
exponential localization of the Wannier functions. As a 
corollary, Wannier functions in Chern insulators are not 
exponentially localized. Moreover, we presented a sim- 
ple criterion to determine when Wannier functions can 
be chosen as real. 

Electron localization is the key of several physical con- 
cepts such as electric polarization js^, piezoelectricity, 
orbital magnetization 37| and the nature of the insulat- 
ing state pj, y] . Our condition for the occurrence of local- 
ization (the vanishing of Chern numbers) is consequently 
a fundamental result for all these subjects. 

We thank K. K. Uhlenbeck, M. Putinar, R. Zentner, D. 
S. Freed, F. Mauri and D. Ceresoli for useful comments, 
and IDRIS (project 061202) for computer time. 
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